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Final  Report  on  the  AFOSR  grant 
No.  FA9550-04-1-0017 

This  Report  presents  a  summary  of  research  conducted  under  the  AFOSR  grant 
No.  FA9550-04-I-0017  “ Atomistic  Modeling  of  Advanced  Intermetallic  Alloys"  (activity  pe¬ 
riod  1 2/2003- 1 1  /2006)  The  goal  of  the  project  was  to  perform  a  fundamental  study  of  structure, 
thermodynamic  and  kinetic  properties  of  structural  intermetallic  alloys  by  means  of  atomistic 
computer  simulations.  W  alloys  of  the  Ni-Al  system  were  chosen  as  model  materials.  The 
Ni3Al-based  7'  phase  is  a  promising  structural  material  for  high-temperature  applications.  Two- 
phase  7-7'  alloys  serve  as  prototypes  of  Ni-base  superalloys.  The  methodology  used  in  this 
project  includes  the  construction  of  new  interatomic  potentials,  large-scale  molecular  dynamics 
(MD)  and  Monte  Carlo  simulations  and  other  methods.  The  properties  studied  here  include 
diffusion,  structure  and  energy  of  generalized  stacking  faults  (GSFs)  in  the  7'  phase,  7/7'  inter¬ 
phase  boundaries,  and  dislocation  behavior  in  the  7'  phase. 

1  Interatomic  potentials  for  the  Ni-Al  system 

A  new  interatomic  potential  has  been  constructed  within  the  embedded-atom  method  (EAM) 
formalism  by  fitting  to  both  experimental  and  first-principles  data  [1].  The  potential  set  includes 
accurate  EAM  potentials  for  Ni  and  A1  and  a  cross-potential  fit  to  first-principles  energies  of 
various  Ni-Al  compounds,  including  L 1 2-Nia A I.  The  potentials  have  been  tested  by  calculating 
a  large  variety  of  thermodynamic  and  physical  properties  of  Ni,  A1  and  Ni3AI,  including  phonon 
frequencies,  thermal  expansion,  melting  temperatures  and  heats,  planar  fault  energies,  point 
defect  formation  and  migration  energies,  as  well  as  experimental  and  first-principles  formation 
energies  of  a  series  of  Ni-Ai  compounds.  The  potentials  has  shown  excellent  performance 
for  most  properties  tested.  The  Ni-rich  part  of  the  Ni-Al  phase  diagram  has  been  calculated 
by  the  thermodynamic  integration  method  in  conjunction  with  grand-canonical  Monte  Carlo 
simulations.  The  reasonable  agreement  with  experiment  (Fig.  1)  makes  this  potential  applicable 
to  atomistic  simulations  of  7-7'  alloys. 

2  Diffusion  in  NiAI  and  Ni3Al 

Diffusion  mechanisms  operating  in  Ni3AI  and  NiAI  have  been  established  by  large-scale  MD 
simulations  [1,2].  Ni  atoms  diffuse  in  Ni3AI  by  exchanges  with  vacancies  on  the  Ni-sublattice, 
whereas  A1  goes  to  the  Ni  sublattice,  forms  Atjy;  antisites,  and  diffuses  by  exchanges  with  Ni 
vacancies.  The  diffusion  mechanisms  in  NiAI  are  more  complex.  In  this  compound,  devia¬ 
tions  from  stoichiometry  on  the  Al-rich  side  are  accommodated  by  structural  vacancies  on  the 
Ni  sublattice,  whereas  the  Ni-rich  side  is  dominated  by  NUt  antisites  [3].  It  was  found  that 
the  atomic  configuration  arising  after  a  nearest-neighbor  jump  of  a  Ni  vacancy  is  mechanically 
unstable.  Because  of  this  instability,  the  vacancy  implements  two  sequential  nearest-neighbor 
jumps  as  one  collective  transition  involving  two  atoms.  Such  collective  jumps  initiate  and  com¬ 
plete  six-jump  vacancy  cycles  of  a  Ni  vacancy.  Thus,  such  cycles  occur  by  either  four  or  three 
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Figure  1:  Ni-Al  phase  diagram  calculated  with  the  EAM  potential  developed  in  this  project 
(points)  in  comparison  with  experiment  (lines). 


collective  vacancy  jumps,  depending  on  whether  the  central  configuration  of  the  cycle  is  me¬ 
chanically  stable  (as  for  [100]  cycles)  or  unstable  (as  for  [110]  cycles).  Calculations  show  that 
[110]  cycles  are  energetically  more  favorable  and  should  happen  more  frequently  than  [100] 
cycles.  Next-nearest-neighbor  vacancy  jumps  have  diffusion  rates  comparable  to  experiment, 
suggesting  that  Ni  diffusion  on  its  own  sublattice  is  an  important  diffusion  mechanism  in  NiAl. 

These  findings  have  been  verified  by  both  static  calculations  and  MD  simulations  [2].  Diffu¬ 
sion  cocffients  ofNi  and  Al  in  Ni:jAI  and  NiAl  by  different  mechanisms  have  been  calculated  by 
the  kinetic  Monte  Carlo  method  as  well  as  by  direct  MD  simulations.  The  calculated  diffusion 
coefficients  compare  well  with  experimental  data.  Together  with  our  previous  AFOSR  work 
of  diffusion  in  Ti-AI  compounds,  we  now  have  a  well  established  methodology  for  predictive 
diffusion  calculations  in  transition-metal  aluminides. 

3  Interfaces  in  7-7'  alloys 

Interfaces  play  an  important  role  in  mechanical  behavior  of  (he  7'  phase  and  7-7'  alloys  of  the 
Ni-Al  system.  The  GSF  energies  dictate  the  dissociated  structure  of  superdislocations  in  Ni3Al 
and  thus  determine  their  mobility.  The  energy  of  7/7'  interfaces  affects  the  coarsening  kinetics 
of  7'  particles  in  the  7  matrix  and  creep  resistance.  Table  I  summarizes  the  interface  energies  at 
T  =  0  K  computed  with  our  EAM  potential.  They  show  very  reasonable  agreement  with  both 
experimental  data  and  first-principles  calculations  [I  j. 

It  was  established  by  grand-canonical  Monte  Carlo  simulations  that  the  (111)  APB  in  stoi¬ 
chiometric  Ni3AI  shows  a  noticeable  antisite  disorder  at  high  temperatures,  whereas  other  GSFs 
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remain  atomically  sharp.  At  high  temperatures,  7/7'  interfaces  are  diffuse,  with  the  long-range 
order  parameter  changing  gradually  from  almost  !  to  zero  over  a  10-15  A  distance  (Fig.  2).  We 
are  currently  in  the  process  of  thermodynamic  calculations  of  interface  free  energies  as  func¬ 
tions  of  temperature  and  chemical  composition  using  thermodynamic  integration  methods  in 
conjunction  with  Monte  Carlo  simulations. 


Table  I :  Energies  (in  mJ/m2)of  generalized  stacking  faults  and  7/7'  interfaces  with  different  orientations 
in  N^Al  calculated  with  the  EAM  potential  [  1 J.  APB:  anti-phase  boundary;  CSF:  complex  stacking  fault; 
SISF:  superlattice  intrinsic  stacking  fault. 
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Figure  2:  Concentration  profile  across  a  (1 00)  7/7'  interface  in  the  binary  Ni- A I  system  obtained 
by  Monte  Carlo  simulations  at  700  K.  Individual  points  represent  the  A1  concentration  averaged 
over  individual  (200)  layers.  The  arrow  marks  the  initial  ( T  =  0)  position  of  the  interface. 


4  The  “prewetting”  transition  in  7-7'  alloys 

Grand-canonical  Monte  Carlo  simulations  of  APBs  in  Ni-rich  y'  reveal  5-10  A  wide  regions 
around  the  APB  with  greatly  reduced  long-range  order  and  chemical  composition  approaching 
the  composition  of  the  7-phase.  As  the  bulk  composition  approaches  the  two-phase  equilibrium 
line,  the  thickness  of  the  7-layer  at  the  APB  rapidly  increases  and  eventually  diverges.  This 
local  disordering  of  the  APB  must  occur  as  a  first-order  interface  phase  transformation  which 
can  be  classified  among  so-called  “wetting”  transitions,  meaning  that  the  APB  becomes  “wet” 
with  the  7  phase.  A  fact  which  is  important  for  applications  is  that  this  transition  can  happen 
within  the  7' -domain  on  the  phase  diagram,  i.e.  before  the  two-phase  equilibrium  is  reached 
(“prewetting”  transition). 
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We  have  studied  the  (1 1 1)  APB  prewetting  transformation  as  a  function  of  temperature  and 
bulk  composition.  Only  Ni-rich  off- stoichiometries  have  been  examined  due  to  their  relevance 
to  7-7'  alloys.  Ni  is  found  to  always  segregate  to  this  APB.  Figure  3  displays  the  excess  amount 
of  Ni  (relative  to  its  bulk  concentration,  measured  in  the  effective  number  of  monolayers)  as 
a  function  of  the  bulk  composition  at  selected  temperatures.  Observe  that  the  amount  of  Ni 
segregation  increases  with  the  degree  of  off-stoichiometry  but  decreases  with  temperature.  At 
a  constant  temperature,  the  segregation  curve  has  two  stages.  At  small  deviations  from  the  sto¬ 
ichiometry,  the  segregation  rapidly  increases  in  a  linear  manner  and  reaches  a  saturation.  At 
this  stage  the  Ni  concentration  within  the  APB  rapidly  grows  whereas  the  APB  width  remains 
almost  constant  and  equal  to  approximately  2d(m).  By  the  end  of  this  stage,  the  chemical  com¬ 
position  within  the  APB  comes  very  close  to  the  equilibrium  bulk  composition  of  the  7-phase. 
At  the  second  stage,  the  width  of  the  APB  region  increases  whereas  its  chemical  composition 
remains  nearly  constant  and  close  to  the  7-phase  composition.  As  the  y  bulk  composition  ap¬ 
proaches  the  miscibility  gap,  the  APB  width  rapidly  grows  and  diverges.  We  are  in  the  process 
of  collecting  enough  data  for  constructing  the  entire  prewetting  line  on  the  Ni-Al  phase  diagram. 
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Figure  3:  Ni  segregation  at  the  (1 1 1)  APB  in  N^Al  alloys  as  a  function  of  the  bulk  composition 
at  selected  temperatures.  The  left  end  of  each  curve  is  close  to  the  bulk  solubility  limit.  The 
arrows  indicate  the  points  of  the  “prewetting”  transformation  which  turns  the  APB  into  a  7 
layer. 


The  APB  prewetting  transition  can  have  significant  consequences  for  the  precipitation  and 
coarsening  of  7'  particles  in  7-7'  alloys.  If  two  growing  particles  can  merge  without  forming  an 
APB,  they  will  readily  do  so.  However,  if  an  APB  should  result  from  the  merger,  such  particles 
will  not  merge  and  will  stay  clear  of  each  other.  This  effect  imposes  an  additional  resistance 
to  the  coarsening  and  is  considered  to  be  responsible  for  the  discontinuous  coarsening/rafting 
effect  and  other  structural  features  of  7-7'  alloys  that  affect  mechanical  properties.  Dislocations 
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in  the  7'  phase  can  also  be  influenced  by  the  APB  prewetting  effect.  Screw  superdislocations 
in  7'  typically  dissociate  on  a  {111}  plane  with  the  formation  of  an  APB  bounded  by  two 
superpartials.  The  APB  replacement  by  a  7  layer  must  have  a  strongly  impact  on  the  dislocation 
mobility  at  high  temperatures.  This  must  affect  the  creep  resistance  of  the  7'  phase  as  well  as 
superalloys. 


5  Dislocations  in  the  7'  phase 

NiaAl  is  one  of  several  Ll2  intermctallic  compounds  which  exhibit  the  yield-stress  anomaly 
(YSA,  the  yield  stress  increases  with  temperature  and  shows  a  peak  at  a  few  hundred  K).  The 
YSA  in  Ni3AI  is  accompanied  by  other  unusual  features  of  plastic  deformation,  such  as  (i) 
strong  orientation  dependence  of  the  CRSS,  (ii)  tension-compression  asymmetry,  and  (iii)  small 
strain-rate  sensitivity.  The  mechanism  responsible  for  this  behavior  is  believed  to  be  thermally 
activated  formation  of  dislocation  locks  arising  by  cross-slip  of  screw  dislocations,  or  their 
segments,  from  the  primary  {III}  plane  to  a  cube  {001}  plane.  After  30  years  of  studies 
of  the  YSA  in  Ni3AI,  there  is  still  a  number  of  unresolved  problems.  For  example,  there  is  no 
consensus  about  the  exact  structure  of  the  dislocation  locks.  The  mechanism  of  the  small  strain- 
rate  sensitivity  is  not  clear  either.  While  dislocations  in  Ni3AI  have  been  studied  by  atomistic 
simulations  in  the  past  [4],  all  previous  work  was  performed  by  static  relaxation  at  0  K.  It  thus 
dealt  mostly  with  postulated  dislocation  structures  and  did  not  directly  address  the  temperature 
effects,  which  are  at  the  heart  of  the  YSA. 

•  m  l  i — * 

(a) 


(b) 


\ 


(0 


(d)  th)  X 

Figure  4:  Superdislocation  configurations  observed  by  MD  simulations.  Thick  lines  -  (111) 
APB,  double  thin  lines  -  (001}  APB,  single  thin  lines  -  CSF,  dots  -  Shockley  partials.  (a) 
Planar  core;  (b)  Nonplanar  core  on  octahedral  planes,  (c),  (d)  PPV  lock  of  one  partial;  (e),  (0 
PPV  locks  of  both  partials;  (g)  Complete  K  W  lock;  (h)  Incomplete  K  W  lock  found  in  this  work. 


We  have  started  an  extensive  study  of  dislocation  behavior  in  Ni3AI  by  means  of  MD  simu 
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lations  using  the  RAM  potential  developed  in  this  project  []].  The  advantage  of  MD  is  that  the 
dislocation  of  interest  evolves  naturally  under  an  applied  shear  stress  and  a  chosen  temperature. 
No  a  priori  assumptions  are  made  about  its  behavior;  instead,  the  mechanisms  of  its  motion 
are  deduced  from  the  simulations.  Our  simulations  use  a  block  containing  up  to  half  a  million 
atoms.  Taking  advantage  of  massive  parallel  computing,  MD  simulation  times  up  to  several  tens 
of  nanoseconds  have  been  achieved.  Several  temperatures  (700-1300  K),  stress  levels  (500-800 
MPa)  and  strain  rates  have  been  tested.  The  initial  configuration  is  always  a  b  =  [llO]  screw 
superdislocation  dissociated  on  the  (111)  plane.  The  shear  stress  is  applied  parallel  to  this  plane 
in  the  [llO]  direction.  The  dislocation  core  structure  and  its  motion  are  constantly  monitored  by 
the  Nye  tensor  method  developed  earlier  in  this  work  (see  below).  The  results  obtained  so  far 
can  be  briefly  summarized  as  follows: 

The  initial  state  of  the  dislocation  is  a  planar  core  dissociated  on  (111)  with  the  partials 
separated  by  an  APB  (Fig.  4).  Each  partial,  in  turn,  is  narrowly  dissociated  in  two  Shockley 
partials  separated  by  a  CSF.  This  configuration  turns  out  to  be  unstable  and  at  finite  temper¬ 
atures  quickly  transforms  to  a  nonplanar  core  with  the  partials  dissociated  on  the  secondary 
octahedral  plane  (111).  This  dissociation  is  dynamic:  every  once  in  a  while  the  Shockey  par¬ 
tials  spontaneously  recombine  and  immediately  redissociate.  Such  recombination-dissociation 
events  require  a  constriction  of  the  CSF  and  are  thermally  activated  events.  Eventually,  one  or 
both  of  the  superpatials  constrict,  cross-slip  along  the  (001)  plane  by  6/2,  and  re-dissociate  on 
a  (111)  plane  parallel  to  the  initial  one,  thus  forming  a  so-called  PPV  lock  [4].  The  dissocia¬ 
tion  plane  of  the  locked  partial  is  again  dynamic  and  can  spontaneously  switch  between  (111) 
and  (111).  This  PPV-locked  configuration  is  more  stable  than  the  unlocked  one,  although  a 
few  spontaneous  unlocking  events  have  been  observed.  The  latter  points  to  a  relatively  weak 
character  of  PPV  locks. 

When  stress  is  applied  to  an  unlocked  dislocation,  it  begins  to  glide  along  the  (111)  plane. 
During  the  “free  flight”  of  the  dislocation  the  separation  of  the  supcrpartials  and  the  width  of 
the  CSFs  are  significantly  smaller  than  in  the  standstill  state  without  a  stress.  After  traveling 
a  distance  of  a  few  tens  of  nanometers,  the  trailing  partial  forms  a  PPV  lock  and  stops  the 
dislocation.  The  dislocation  remains  locked  for  a  few  hundred  picoseconds  before  the  trailing 
partial  breaks  the  lock  and  the  glide  continues.  Breaking  a  PPV  lock  requires  a  constriction 
of  the  Shockley  partials  and  a  cross-slip  by  6/2  along  the  (001)  plane.  Overall,  the  disloca¬ 
tion  moves  in  a  stop-and-go  manner  (“jerkey  motion”)  caused  by  repeated  PPV  locking  and 
unlocking  events  (Fig.  5).  Although  the  average  dislocation  velocity  (100-200  m/s)  is  much 
smaller  than  for  the  velocity  of  free  glide,  the  resistance  imposed  by  the  PPV  locks  is  not  strong 
enough  to  account  for  the  YSA.  Furthermore,  this  resistance  decreases  with  temperature  since 
the  unlocking  process  is  assisted  by  temperature. 

After  a  few  nanoseconds  of  stop-and-go  motion,  the  dislocation  stops  due  to  the  formation  of 
a  much  stronger  lock  than  PPV,  which  blocks  its  motion  for  the  rest  of  the  simulation.  Namely, 
in  the  process  of  unlocking  of  a  PPV  configuration,  the  constricted  trailing  partial  “misses” 
the  primary  (111)  plane  and  continues  to  glide  along  the  (001)  cube  plane  by  a  few  more 
Burgers  vectors  before  it  stops  and  dissociates  into  Shockley  partials  in  an  octahedral  plane. 
Subsequently,  it  may  advance  along  the  (001)  plane  somewhat  further  whereas  the  leading 
partial  continues  to  be  on  the  (111)  plane.  This  creates  a  highly  stable  angular  configuration. 
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Figure  5:  Position  of  the  (110)  screw  superdislocation  as  a  function  of  time  in  MD  simulation 
at  the  temperature  of  1000  K  under  a  shear  stress  of  700  MPa.  The  dislocation  first  moves  by 
a  stop  and  go  mechanisms  by  the  formation  and  destruction  of  PPV  locks.  After  1 .5  ns  it  stops 
due  to  the  formation  of  an  incomplete  KW  lock  and  remains  immobilized  until  the  end  of  the 
simulation. 


It  is  pulled  forward  by  the  Peach-Koehler  forces  acting  on  both  partials  but  is  pinned  by  the 
intersection  line  of  the  (11 1)  and  (001)  APBs.  This  configuration  can  only  unlock  if  the  leading 
partial  moves  back  to  the  APB  intersection,  which  requires  work  against  the  Peach-Koehler 
force. 

This  sessile  configuration  can  be  viewed  as  a  new  form  of  a  Kear-Wilsdorf  (KW)  lock.  Its 
formation  is  associated  with  a  large  thermal  activation  (multiple  constrictions  of  the  CSF).  Our 
simulations  suggest  that  this  lock  may  be  responsible  for  the  immobilization  of  dislocations 
in  Ni3Al  and  thus  the  increase  of  its  strength  with  temperature.  The  loss  of  strength  beyond 
the  peak  temperature  is  likely  to  be  caused  by  the  thermally  activated  unlocking  of  such  con¬ 
figurations  and  initiation  of  cubic  slip,  In  the  ongoing  work,  we  are  continuing  to  study  the 
dislocation  locking  mechanisms  in  a  more  systematic  manner  by  analyzing  the  trends  across 
different  temperatures  and  stresses. 

6  The  Nye  tensor  method  for  dislocation  characterization 

A  new  method  of  presentation  of  dislocation  core  structure  by  a  distribution  of  its  Nye  tensor 
components  has  been  developed  [5],  This  method  is  based  on  the  concept  of  a  continuously 
dislocated  crystal  which  represents  the  misfit  due  to  a  dislocation  by  a  distribution  of  infini¬ 
tesimal  dislocations  [6-8].  The  infinitesimal  Burgers  vector  db  of  the  continuously  distributed 
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dislocations  crossing  a  small  area  dS  with  a  unit  normal  n  is  given  by  dh  =  a  ■  n  dS,  a  being 
the  Nye  tensor.  A  computer  program  has  been  developed  that  calculates  the  Nye  tensor  dis¬ 
tribution  over  a  chosen  plane  given  atomic  coordinates  in  a  dislocated  crystal  and  the  perfect 
lattice  orientation.  This  distribution  is  represented  by  contour  plots  of  individual  Nye  tensor 
components.  The  peaks  of  the  distribution  represent  full  or  partial  dislocation.  The  sign  of 
the  relevant  Burgers  vector  component  can  be  immediately  determined  from  the  sign  of  the  re¬ 
spective  component  of  a.  Furthermore,  by  integrating  the  distribution  of  a  over  a  peak  area  all 
components  of  the  Burgers  vector  of  the  dislocation  can  be  determined  with  high  accuracy.  It 
has  been  demonstrated  that  this  method  is  insensitive  to  lattice  distortions  other  than  those  as¬ 
sociated  with  dislocation  cores,  which  makes  its  spatial  resolution  much  higher  than  with  other 
methods.  As  an  illustration,  Figure  6  shows  the  structure  of  a  {110)  screw  superdislocation  in 
Ni3AI  dissociated  with  the  formation  a  (111)  APB.  Each  superpartial  is  observed  to  dissoci¬ 
ate  into  closely  separated  Shockley  partials  on  a  non-coplanar  {111}  plane.  The  plot  does  not 
display  the  GSFs  in  order  to  enhance  the  clarity  of  visualization  of  the  partial  dislocations. 

The  Nye  tensor  method  has  been  applied  to  the  intepretation  of  high-resolution  electron 
microscopy  images  of  screw  dislocations  in  Mo  [9].  It  has  been  demonstrated  that  Nye  tensor 
plots  of  the  images  “wash  out”  the  effect  of  the  Eshelby  twist  and  give  a  more  informative 
presentation  of  the  core  structure  than  other  methods  can  do. 
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Figure  6:  Nye  plot  of  a  (110)  screw  dislocation  in  Ni3AI  after  a  molecular  dynamics  run  at  500 
K.  The  contour  plot  shows  the  Nye  tensor  component  parallel  to  the  dislocation  line  and  to  the 
[112]  direction.  The  dashed  line  was  drawn  by  hand  to  outline  the  (111)  APB.  The  brown  and 
blue  peaks  identify  the  Shockley  partials  and  indicate  the  signs  of  their  [112]  components.  The 
distances  are  in  A  and  the  Nye  tensor  in  1/  A. 
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